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Galactic bulge formation as a mctximum intensity starburst 

Bruce G. Elmegreen^ 

ABSTRACT 

Properties of normal galactic star formation, including the density dependence, threshold 
density, turbulent scaling relations, and clustering properties, are applied to the formation of 
^ , galactic bulges. One important difference is that the bulge potential well is too deep to have 

I allowed self-regulation or blow-out by the pressures from young stars, unlike galactic disks or 

dwarf galaxies. As a result, bulge formation should have been at the maximum rate, which is 
such that most of the gas would get converted into stars in only a few dynamical time scales, 
or ~ 10^ years. The gas accretion phase can be longer than this, but once the critical density 
^ ■ is reached, which depends primarily on the total virial density from dark matter, the formation 

•/^ , of stars in the bulge should have been extremely rapid. Such three-dimensional star formation 

should also have formed many clusters, like normal disk star formation today. Some of these 
clusters may have survived as old globulars, but most got dispersed, although they might still 
, be observable as concentrated streams in phase space. 

:^ 

o 

^ ' 1. Introduction 

^ ' On a large scale, star formation in galaxies appears to be regulated by the self-gravity of the gas, 

producing new stars at a rate proportional to the average gas density p multiplied by the standard 

■ ■ 1/2 

, gravitational rate, {Gp) . This leads to the Schmidt law of star formation (Madore 1981; Larson 1988; 

^ I Elmegreen 1991; Kennicutt 1998) and the observation of a decreasing gas consumption time at higher 

densities, which is essentially the starburst phenomenon (Larson 1987; Elmegreen 1994; Kennicutt 1998). 

Large-scale star formation also has a density threshold (Madore, van den Bergh, & Rogstad 1974; 
Gallagher & Hunter 1984; Guiderdoni 1987), below which the rate per unit gas mass is much less than 
the inverse gravity time. This threshold has been attributed to several types of limits: (1) a minimum 
column density for strong self-gravity of a gas disk in the presence of galactic differential rotation, 
i.e., acrit ^ Ka/ {ttG) (Toomre 1964; Goldreich & Lynden-Bell 1965; Quirk 1972; Fall & Efstathiou 
1980; Zasov & Simakov 1988; Kennicutt 1989; Elmegreen 1991), and (2) a minimum column density 
corresponding to the minimum pressure (which scales as ^ Gcr^) that gives a cool thermal gas phase, 
iycrit,th = 7(tTg/ [Oar])^'^ (///q)"'^ M© pc^^ for radiation field / relative to the local value, /©, total 
gas-|-stellar column density ctt, and ratio of thermal to total gas pressure, 9 (Elmegreen & Parravano 1994). 
The outer edges of galaxies where star formation slows abruptly (Ferguson et al. 1998) are probably the 
result of contributions from both factors, considering that gaseous self- gravity becomes low (Kennicutt 
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1989) and the HI becomes nearly pure warm- phase (Braun 1997). In addition, the outer parts of galaxies 
are beyond the realm of strong stellar density waves, particularly outside the outer Lindblad resonance for 
the dominant pattern speed of spiral wave modes in the disk (Bertin et al. 1989). The absence of stellar 
density waves removes one of the primary mechanisms for organizing the gas into star- forming clouds. 

On smaller scales, there is a correspondence between the duration of star formation in a region and 
the turbulent crossing time (Efremov & Elmegreen 1998), which suggests that self-gravity and turbulence 
together partition the gas into clouds and cloud clumps, and that the clouds are rather transient, turning 
into stars quickly before a chaotic turbulent boundary destroys the favorable environment (e.g., see 
turbulence simulations in Ballesteros-Paredes, Vazquez-Semadeni & Scalo 1998). The observation of similar 
initial stellar mass fimctions in clusters (Massey 1998) and in integrated galactic light (Bresolin & Kennicutt 
1997) is consistent with this turbulent time scale, because it implies that intermediate and high mass 
stars are not the primary agents of destruction for molecular clouds. Such destruction would steepen the 
combined IMF from all clouds above the Salpeter power-law slope of ~ —1.35 to a slope of ~ —4 (Elmegreen 
1998, 1999), which is observed only rarely in remote regions of the field (Massey et al. 1995). 

Also consistent with this short time scale is the common observation that yomig or embedded star 
clusters still contain significant subclustering, presumably from the pre-stellar cloud. All of this stellar 
substructure would have been mixed together if the cluster age were much more than a crossing time. 
Examples of young clusters with small-scale subclustering are W33 (Beck et al. 1998), G 35.20-1.74 (Persi 
et al. 1997), IC 348 (Lada & Lada 1995), NGC 2264 (Piche 1993), and many other clusters in the Large 
Magellanic Clouds (Elson 1991). 

These observations are useful for understanding the connection between large and small scale star 
formation. Stars form in self- gravitating clouds on a wide range of scales, but the density of these clouds 
increases, and the gravitational time scale decreases, for smaller scales, probably because of the way 
turbulence structures the gas. This structure is also hierarchical, so many small regions of star formation 
come and go while the large region enclosing them continues to form stars. As a result, the large scale 
controls the overall star formation time, while the smaller scales determine the net efficiency, i.e., the 
fraction of the gas that is turned into stars in all of the transient small scale events. 

The surface density threshold for star formation has been observed for a wide variety of regions, 

including normal spiral galaxies (Kennicutt 1989; Caldwell et al. 1992), low surface brightness galaxies (van 
der Hulst et al. 1993), starbursts (Shlosman & Begelman 1989; Kennicutt 1998) and embedded gas disks 
in elliptical galaxies (Vader & Vigroux 1991). However, in all cases considered so far, the geometry has 
been either disk-like, as for large parts of galaxies, or ring-like, as in some nuclear starbursts (Elmegreen 
1994). A threshold for star formation in spheroidal systems like dwarf galaxies and galactic bulges has 
not been seriously considered. Nevertheless, the basic principles of star formation that have been derived 
for disk galaxies should also have applications in other regions, including the galactic bulge, as long as 
the three-dimensional geometry for self-gravity and various turbulent and magnetic processes are properly 
considered. 

Here we discuss theoretical and observational evidence that the criterion for star formation underlying 
most of the galactic observations is one in which the volume density of gas is comparable to or larger than 
the total virial density in that region, including stars and dark matter. This threshold has already been 
suggested for dwarf and irregular galaxies (Elmegreen et al. 1996). It also has a natural extension to three 
dimensions, and therefore should be appropriate to the formation of galactic bulges. In addition, we suggest 
that the likely formation time of the bulge was only a few dynamical times, as it is in local molecular 
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clouds, and that this short time leads to bulge formation luminosities that are extremely large, consistent 
with observations of galaxies at high z (Giavalisco et al. 1995, 1996; Steidel et al. 1996). Moreover, the 
pressures resulting from star formation give an interesting condition for self-regulation and blow-out that 
depends only on the depth of the potential well. This condition is such that normal galaxy disks and dwarf 
galaxies can regulate their star formation to remain near or below the threshold, but galactic bulges and 
whole elliptical galaxies cannot. 

Another implication of using normal star formation processes for galactic bulges is the prediction that 
stars should form in numerous dense clusters. The remnants of these clusters might be visible today as 
stellar streams or other concentrations in phase space, and as bulge (and halo) globulars. 



Spitzcr (1942) noted the existence of a maximum stable gas mass for a given temperature in a fixed 
background potential from stars. This maximum mass corresponds to a threshold ratio of the gas density 
to the virial density. Above this density, gas pressure cannot support a cloud against the combined gravity 
from itself and the background stars. 

Some applications of this three-dimensional threshold density have already been made for dwarf and 
irregular galaxies. The sizes of the largest regions of star formation ("star complexes") in dwarf, irregular, 
and spiral galaxies spanning a factor of 10^ in luminosity were found to always equal half the Jeans' length 
for a fixed gaseous velocity dispersion of ~ 10 km s~^, provided the gas density is comparable to the 
background virial density that is given by the TuUy-Fisher relations between galaxy size, dynamical velocity, 
and absolute magnitude (Elmegrccn et al. 1996). This implies that when the gas density exceeds the virial 
density, star formation begins on a scale equal to the Jeans length. These results are interesting here 
because dwarf galaxies are not thin disks and they have no spiral density waves. Thus the star formation 
threshold for dwarfs might be similar to that for other three-dimensional structures. 

Spitzer's derivation was based on the virial theorem. It gave the result that an isothermal gas cloud 
embedded in a stellar spheroid is unstable if its central density Pc exceeds Pvir(l + /3/2)~^ for gaseous 
gravitational potential pGMgas/ R- 

We recalculate this result here by direct integration of the equations of hydrostatic equilibrium 
for a mixture of cool gas and a dissipationless fluid that might be considered to be stars, dark matter, 
or both. Here we just call the second component dark matter. The equations of equilibrium are 
dPgas/dr = —pgas9 and dPjjM/dr = —p^Mg, respectively. The gravitational acceleration comes from 
both components according to Poisson's equation in spherical coordinates: V • g = 47rG {pgas + Pdm)- 
The results are normalized to the central dark matter density, poM o, the dark matter velocity dispersion, 

1/2 ' 

dDM = (Pdm/ Pdm) , and the gravitational length and mass scales given by these quantities, namely. 



The top left diagram of figure 1 plots various density profiles for dark matter (dashed lines) and gas 
(solid lines), considering different central gas densities. The lower dashed line corresponds to the upper 
solid line. The top right diagram shows the total gas mass as a function of the central gas density, for 
various ratios of the gaseous velocity dispersion to the dark matter dispersion. When the gaseous velocity 

dispersion is small, the total gas mass peaks at intermediate values of the central gas density. This is 
because at lower gas density the gas sphere has a low density throughout a large size, and at higher gas 
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Pz3M,o)' respectively. 
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density, the gas sphere has a high density throughout a small size. As the gaseous dispersion approaches 
the stellar dispersion, however, the total mass of the gas sphere approaches infinity, as does the mass of the 
stellar sphere, because of the usual infinite extent of isothermal sphere solutions. 

The lower left diagram in figure 1 shows the peak total gas mass versus the velocity dispersion ratio 
as a solid line, and the ratio of the total gas mass to the dark matter mass inside the gas sphere as a 
dashed line. The gas mass ratio is about 0.03 at low dispersion ratio. The final result is the gas density 
ratio, plotted in the lower right. When the central gas density exceeds about 12.6 times the central dark 
mass density, or the average gas density inside the half-mass gas radius exceeds about 6 times the central 
dark matter density, there are no stable solutions which permit the addition of more mass to a cool gas 
component of a gas+dark matter galaxy. Any additional gas mass leads to catastrophic collapse of the 
whole gas sphere. If write the total gas+dark matter density as the virial density, then this result gives a 
critical average gas density 

Pcrit Pvir- (1) 

This is about the same value as the critical gas density in spiral disks. The disk condition, from 
Kennicutt (1989), is usually written in terms of the surface density, 

'-"-"■'St- 

For a flat rotation curve with speed V and radius R, the epicyclic frequency k = To get the virial 

density, wc take V"^ = GM/R = AiiGpyirR^ /'i- The surface density gets converted into a volume density 
using the scale height H, as p = a {2H)^^. The scale height for gas is about H = Ugas {2TrGpT)~^^^ , where 
Pt is the total midplane density, including stars. Thus Pcrit,disk ^ 0.8 (pvirPr)^^^ ■ Locally, p^i^ ~ 2.5 x 10~^^ 
g cm~^ for V = 220 km s^^ and R = 8.5 kpc; this is about 1/3 the total midplane density, so we take 
Pt ^ ^Pvir as representative. Then 

Pcrit,disk ^ ^-^Pviv (3) 

The virial density is a reasonable threshold for major events of galactic star formation because the 
main processes that resist gaseous self-gravity in disks, namely, shear and rotation, have effective time scales 
(Oort and respectively) that equal approximately {Gpvir) Gaseous self-gravity, on the other 
hand, has an effective time scale of {Gpgas)~^^^ ■ Thus self-gravity dominates when pgas > Pvir over a length 
scale larger than the Jeans length. The interstellar magnetic field will not change this result much if it is in 
equipartition with the turbulent motions, because the first step in cloud formation following a gravitational 
instability consists of motions along the magnetic field, and the field pressure does not resist this motion 
much. Spiral waves in disks may trigger substantial amounts of star formation, even in sub-critical disks, 
but since the wave incidence rate is on the order of a rotation rate, the net time scale for star formation, 
averaged over a rotation, is still about the orbit time (Wyse & Silk 1989). In addition, spiral waves seem to 
require dense gas for strong amplification (Lubow, Cowie, & Balbus 1986), so again we get the result that 
when the gas density is near the threshold for strong self-gravity, at Pgas Pvir^ star formation, even in 
spiral density waves, becomes rapid. 



3. The threshold stctr formation rate 



Considering the above discussion, the threshold star formation rate per unit volume in any region will 
be about equal to the threshold gas density, which is essentially the local virial density, divided by the local 
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dynamical time, or 

SFR^f^^eG'/-'pT^^' (4) 

'-dyn '^^dyn 

where e 0.1 is the fraction of the gas that is turned into stars in each major star formation event (given 
by the combined efficiencies of star formation on smaller scales). The dynamical time is tdyn ^ R/V. For 
tdyn — lO^tdyn,? years, this rate becomes 

SFR « 4^ kpc"^ yr-^ (5) 

^dyn,7 

Meurer et al. (1997) found an upper limit to the areal star formation rate in starburst nuclear disks 
that is equal to ~ 45 Mq kpc~^ yi'~^. This comes from equation ^ if we multiply the volume star formation 
rate by a galactic thickness of ^ 0.3 kpc, use e ~ 0.1, and take a dynamical time appropriate for a nuclear 
disk, tdyn = R/V ^ 5 My. Equation || predicts that even higher star formation rates will be found in 
smaller regions where tdyn is smaller, although the fraction of galaxies undergoing such intense bursts will 
be small because of rapid gas consumption during the active phase. 

Meurer et al. (1997) also discussed wind limitations to star formation rates. They estimated the 
pressure from a star-forming region to be 

PSF 7.2 X 10-"dyne cm-^ x ( ^ (6) 

VM0kpc-V-V 

= 1.1 X lO^ne cm-2 x ( . (7) 

\gmcm ^s / 

This pressure is proportional to the two-dimensional star formation rate, SFR2D, which is the galactic 
thickness 2H times the 3-dimensional rate from equation ^. Thus, at threshold (in cgs units), 

PsF ~ 1.1 X lO^e [Gp^^rf^ Pv^r2H. (8) 

This disruptive pressure should be compared with the self-gravitational pressure that holds 
the region together, which is proportional to the square of the mass column density. For a disk, 
Pgrav,disk ~ (7''/2) GagasO'T for total column density in the gas layer, ar- This equation follows from the use 
of = flgjjg/ {2ttGpt) for gas disk scale height H and from the relation <j = 2Hp; then Pgrav = PgasO^aai 
as expected. If we take pgas ~ Pvir at threshold, and px ^ Spvir as above, then Pgrav, disk ~ GnGpl^j.H'^ . 
For a sphere, Pgrav, sph (ti'/Z) GpTPgas'^H^ again if we define H = agas/ i^TrGpr)^^^ in the same way. But 
now pt = Pvir for the problem considered in figure 1, and Pgas ^ Pvir inside the half-mass radius. Thus 
Pgrav, sph ~ 27rG'/9^j^-ff ^ . To covcr both cases, we write this pressure as 

Pgrav = anGpl^^H^, (9) 

where a = 6 and pT = ipvir for a disk, and a — 2 and pT — pvir for a sphere. 

Now we see that the star formation pressure exceeds the self-gravitational binding pressure when 

1.1 X lO^e > (7r/2) a {Gp^,„f'^ H ^ llagas] (10) 

the coefficient 1.7 is the average between (I.Stt)^^^ — 2.2 for a disk and (O.Stt)^^'^ = 1.2 for a sphere. The 
coefficient on the left is in cgs units, namely, cm s~^. Thus star formation pressure can blow apart a region 
at the threshold density if the gaseous velocity dispersion, in km s^^, is less than a certain limit, 

o,gas < 65e km^^ s^^. (11) 
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This limit is surprisingly simple: for typical e 0.1, it says that star formation can regulate or stop 
itself in a disk (e.g., Franco & Shore 1984) or dwarf galaxy (e.g., Dekel & Silk 1986) where the gaseous 
velocity dispersion is low, Ugas 7 km s~^, but it cannot regulate or stop itself in a young galactic bulge 
or elliptical galaxy, where the velocity dispersion from the potential well is large, Ugas > 100 km s^^. Star 
formation might also fail to regulate itself in nuclear disks or interacting galaxies, which have large gaseous 
velocity dispersions too (Elmegreen et al. 1993; Irwin 1994; Elmegreen et al. 1995; Brinks et al. 1997; 
Kaufman et al. 1997). Such failure to self-regulate should lead to rapid and intense star formation on a 
time scale of only a few dynamical times. 

Rapid star formation in bulges and elliptical systems also leads to large luminosities. We convert the 
star formation rate from equation [s] to a volume emissivity jsF, following Meurer et al. (1997), and get 



9 X 10"e , , _3 
JSF 73 L© kpc . (12) 

''dynj 



This expression may be used to get the total luminosity by multiplying it by the volume, 47ri?'^/3 for radius 
R, and then the result can be written as a function of the dynamical speed, V, using tdyn = R/V . The 
resultant luminosity is 

^-4^io"KiooiL-t)'^- 

Evidently, the threshold starburst luminosity also depends only on the depth of the potential well of 
the galaxy. This result should apply to many types of systems, since we have assumed very little about the 
nature of the star formation or how the gas got into the potential well. For merging galaxies, the dynamical 
speed can be '--^ 200 km s~^ or more, leading to star- formation luminosities of L 10^^ for times of 
idyn ^ 10* years {tdyn/^ ^ 10^ years if the residual gas from inefficient star formation continues to make 
stars). For bulges in which V ^ 100 km s^^, the luminosity during bulge formation can be L 10^^ Lq for 
a total formation time of tdyn/^ ^ lOR/V ^ 10* years. These values are consistent with observations of the 
Lyman a emitting galaxy at z — 3.4 found by Giavalisco et al. (1995), who suggest they are witnessing 
the formation of a galactic bulge or elliptical galaxy (see also Lilly 1999). According to this theory, the 
formation of a typical galactic bulge could be all over within ~ 100 million years once it starts, i.e., within 
~ tdyn/s for a complete conversion of gas into stars. This result is consistent with recent observations of 
rapid bulge formation in the Milky Way (Renzini 1999). 

The presence of a critical density for star formation in three-dimensional systems implies that there 
should have been a delay in the onset of bulge formation during part of the accretion and cooling phase of 
the galactic gas near the center. Star formation had to wait until the accretion built up a gas density in 
the center comparable to ~ lOx the density of dark matter there. After this critical density was reached, 
star formation should have been extremely rapid, producing 10^^ L© for 10* years and then dying off 
(depending on the rate of continued gas accretion). The bulge should not have been able to regulate this 
rapid star formation by supernovae, winds or other young stellar pressures, because the total binding 
pressure exceeded the star formation pressure. Any winds that were generated had to move between the 
dense gas clumps and find their way out of the bulge in streams without completely expelling these clumps 
from the system. 

This starburst model of bulge formation differs from the model by Wada, Habe & Sofue (1995), who 
suggested that a nuclear starburst ejects a giant gas bubble, which forms stars and relaxes into a bulge. It 
also differs from other bulge formation models in which the bulge stars are ejected individually from the disk 
as a result of bar-driven or other resonances (Combes et al. 1990; Pfenniger & Friedli 1991; Friedli & Benz 
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1995) , or ejected in bulk when a bar dissolves (Hasan, Pfenniger, & Norman 1993; Friedli & Benz 1993). In 
the present model, the bulge formed in a three-dimensional distribution (Eggen, Lynden-Bell, and Sandage 
1962) as it does in the clusters of local molecular clouds, except on a much larger scale. A prediction of 
this model can then be made by analogy to all other known star formation: the bulge should have been 
born with an extremely clumpy distribution of young stars, i.e., composed almost entirely of clusters with a 
hierarchical distribution and a mass function something like M~'^dM (Fleck 1996; Elmegreen & Falgarone 

1996) . Presumably some of these clusters remain today as globular clusters (Surdin 1979; Okazaki & Tosa 
1995; Elmegreen & Efremov 1997), but most should have been dissolved by interactions and variable tidal 
forces (Gncdin & Ostrikcr 1997). In that case, they may still be observable as stellar streams, which arc 
concentrations in phase space. This was apparently the fate of local disk star clusters that were dispersed 
with age (Eggen 1989; Dehnen 1998). 
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Fig. 1. — (top left) Density profiles for gas and dark matter in a spherical potential, considering different 
central gas densities, (top right) Gas mass versus central density for four values of the ratio of the gas to 
the dark matter velocity dispersion, (lower left) Maximum gas mass, and ratio of the maximum gas mass to 
the dark matter mass in the gas sphere, versus velocity dispersion ratio, (lower right) Central gas density at 
the maximum stable gas mass (solid line) and average gas density inside the radius containing half the total 
gas mass (dashed line), versus velocity dispersion ratio. 



